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Abstract. Recent results indicate that the Babcock-Leighton mechanism for poloidal 
field creation plays an important role in the solar cycle. However, modelling this mech- 
anism has not always correctly captured the underlying physics. In particular, it has 
been demonstrated that using a spatially distributed near- surface alpha-effect to para- 
metrize the Babcock-Leighton mechanism generates results which do not agree with 
observations. Motivated by this, we are developing a physically more consistent model 
of the solar cycle in which we model poloidal field creation by the emergence and flux 
dispersal of double-rings structures. Here we present preliminary results from this new 
dynamo model. 



1. Introduction 



It is observed that the number of sunspots on the solar surface varies periodically on the solar 
surface with a periodicity of 1 1 years. One of the most promising approaches for modelling 
this solar cycle is using axisymmetric kinematic solar dynamo models. In such models, toroidal 
field generation takes place in the tachocline by stretching of poloidal field lines in <p direction 
by the differential rotation, while poloidal field generation takes place due to decay of tilted 
bipolar sunspot pairs near solar surface, commonly known as Babcock-Leighton mechanism. 
The widely used approach of treating the Babcock-Leighton mechanism by a spatially distributed 
near-surface alpha effect is computationally convenient but has some inherent problems. Here we 
develop a model for the solar cycle, where we model the generation of poloidal fields by the 
emergence and flux dispersal of double ring structures. 
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2. Model and Results 

Here we assume axisymmetry in our calculations. The axisymmetric dynamo equations are 

d A + I [ Vp . V(M)] = 77 (v 2 - 4) A + *o/(r, 0)F{B tc )B tc (1) 
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where A is the poloidal field, B is the toroidal field , v p is the meridional flow, Q is the differential 
rotation, ?] is the turbulent magnetic diffusivity and s = r sin 6. For simulations with double ring 
algorithm, Kq = 0. 

In this model we use a two step radially dependent magnetic diffusivity profile as described in 
Munoz-Jaramillo, Nandy & Martens (2009). In our case, the diffusivity at the bottom of convec- 
tion zone is 10 8 cm 2 /s, diffusivity in the convection zone is 10 1 l cm 2 /s and supergranular diffus- 
ivity is 5 x I0 l2 cm 2 /s. r cz = 0.73fl o , d cz = O.O25fl , r sg = 0.95fl o and d sg = 0.015fl Q describes 
the transition from one diffusivity value to another. We also use same differential rotation profile 
as described in Munoz-Jaramillo, Nandy & Martens (2009). 

We generate the meridional circulation profile (v p ) for a compressible flow inside the convection 
zone from this equation V.(pv p ) = so, pu p = V x (i/^) 
where i// is 

{ 1 - e-l 3 ' r0e }{ 1 - eft^/ 2 ) } e -(("o)/r) 2 (3) 

Here i/>o is the factor which determines the the maximum speed of the flow. We use the following 
parameter values B x = 1.5, B 2 = 1.8,e = 2.0000001, r = (R - R b )/4,T = 3.47 x 10 8 ,y = 
0.95, m = 3/2. Here R p = 0.65/?o is the penetration depth of the meridional flow. In this present 
work we use the surface value of meridional circulation as 14m/sec. 



2.1 Modelling Active Regions as Double Rings 

It is believed that poloidal field generation inside the sun takes place through the Babcock- 
Leighton effect and helical turbulence. Here we only consider the Babcock-Leighton alpha effect, 
wherein, poloidal field generation takes place due to decay of bipolar sunspot region. To model 
active region emergence, we follow the methods of double ring proposed by Durney (1997), 
Nandy and Choudhuri (2001), Munoz-Jaramillo, Nandy, Martens & Yeates (2010) and Nandy, 
Munoz-Jaramillo & Martens (2011). We define the (p component of potential vector A corres- 
ponding to active region as 

A ar {r,6) = KiA(<b)F{r)G(!S), (4) 
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Figure 1. (a) Diagram illustrating the quantities which define the latitudinal dependence of a double-ring bipolar pair, 
(b) Poloidal field lines of our double-rings, one is in northern hemisphere and other is in southern hemisphere. 

where K\ is a constant which ensures super-critical solutions and strength of ring doublet is 
defined by A(<D) . We define F(r) as 



r . , _ I r<R -R ar 

'^[^-(fle-^))] r>R e -R ar ' W 



sin 



where R G is the solar radius and R al = O.85/? corresponds to the penetration depth of the AR. 
We define G(6) in integral form as 

G ^=— a f [B-(0') + B + (O')]sm(e')de', (6) 
sin J,, 

where B + (B ) defines the strength of positive (negative) ring: 



B ± (0) = 



< B ar T f - f 

[l + COS -0 ar ±f))] flf . + f-f <0<0 flr + f + f • (7) 

> e ar + f + f 



Here ar is emergence co-latitude, A is the diameter of each polarity of the double ring and 
X = arcsin[sin(y) sin(A fl) -)] is the latitudinal distance between the centres, where the angular 
distance between polarity centres A ar = 6° and the AR tilt angle is y. We set A, i.e diameter of 
ring doublet as 6°. 



2.2 Recreating the Poloidal Field 



In this model, first we check where the toroidal field is higher than the buoyancy threshold at the 
bottom of convection zone in both northern and southern hemispheres. Then we choose one of the 
latitude randomly from both the hemisphere simultaneously at a certain interval of time, using a 
non uniform probability distribution function such that randomly chosen latitudes remain within 
the observed active latitudes. The probability distribution function is made to drop steadily to zero 
between 30° (-30°) and 40° (-40°) in the northern (southern) hemisphere. Second we calculate 
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the magnetic flux of this toroidal ring. Then we find tilt of corresponding active region, using the 
expression given in Fan, Fisher & McClymont (1994) 

yKQ> x /4 B- 5/ \m(A), (8) 

where Bo is the local field strength, <t>o is the flux associated with the toroidal ring and A is the 
emergence latitude. We set the constant such that tilt angle lies between 3° and 12°. 
Third, we remove a chunk of magnetic field with same angular size as the emerging AR from this 
toroidal ring and calculate the magnetic energy of the new partial toroidal ring. Then we fix the 
value of toroidal field such that the energy of the full toroidal ring filled with new magnetic field 
strength is the same as the magnetic field strength for the partial toroidal ring. Finally, we place 
the ring duplets with these calculated properties at the near-surface layer at the latitudes where 
they erupt. 

When we solve these equations with appropriate boundary condition, we can successfully 
simulate the 11 year solar cycle and reproduce the properties of the observed solar butterfly 
diagram. 
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Figure 2. The background shows weak diffuse radial field on solar surface. Eruption latitudes are denoted by symbols 
black (o) and red (+), indicating underlying negative and positive toroidal field respectively. 

We anticipate that solar dynamo models based on this double ring algorithm would be quite 
useful in simulating and understanding various features of the solar cycle. 
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